hep-ph/0108259 
CERN-TH/2001-234 
TAUP-2678-01 



On Electron-Positron Annihilation into 
Nucleon-Antinucleon Pairs 



John EllisU and Marek KarlinerfJ 



a ) Theory Division, CERN, Geneva, Switzerland 

b ) School of Physics and Astronomy, 
Raymond and Beverly Sackler Faculty of Exact Sciences 
Tel Aviv University, Tel Aviv, Israel 

Abstract 

We discuss the puzzling experimental results on baryon-antibaryon 
production in e + e~ annihilation close to the threshold, in particular 
the fact that a(e + e~ — > ran) <; a(e + e~ — ► pp). We discuss an inter- 
pretation in terms of a two-step process, via an intermediate coher- 
ent isovector state serving as an intermediary between e + e~ and the 
baryon-antibaryon system. We provide evidence that the isovector 
channel dominates both e + e~ — > pions and from NN annihilation at 
rest, and show that the observed ratio of a(e + e~ — > nn)/a(e + e~ — > 
pp) can be understood quantitatively in this picture. 
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Experimental data from the FENICE collaboration [|IJ indicate that 
cr(e + e~ — > nn) is relatively large close to threshold. Their data may be com- 
pared with earlier data on e + e~ — > pp [0, ^ and also with data on the time- 
reversed reaction pp — > e + e~, for which more precise data are available [f|]. 
As seen in Fig. p], the combined data indicate that cr(e + e~ — > nn)/a(e + e~ — > 
pp) ^ 1 when E C m ~ 2 GeV. Averaging over the available data on both the 
direct and time- reversed reactions, which are very consistent, and ignoring 
any possible variation with energy, we find: 

a(e+e- -> pp) n +0 . 16 



o.66±g : i5 (i) 



The fact that this ratio is less than unity requires confirmation, but even equal 
cross sections for e + e~ — > pp and e + e~ — > nn would be quite surprising. 

We recall that the ratio of the cross sections for the corresponding t- 
channel processes ep(n) — > ep(n) should be infinite at zero momentum trans- 
fer, where the form factor simply measures the total proton and neutron 
charges, corresponding to a coherent sum over the electromagnetic charges 
of their constituent quarks. It is also believed that the ratio (JTT) should be 
large at high momentum transfers. In a naive perturbative description of 
e + e~ annihilation into baryons, the virtual time-like photon first makes a 
'primary' qq pair, which is then dressed by two additional quark- ant iquark 
pairs that pop out of the vacuum. This dressing is thought to be a perturba- 
tive QCD process at high momentum transfers, which does not distinguish 
between the u and d quarks, since gluon couplings are flavor-blind. Thus, in 
this conventional perturbative picture, the only difference between the pro- 
duction rates of proton and neutron is through the different electric charges 
of the primary qq pairs. The total perturbative cross section is obtained by 
superposing the amplitudes with different primary qq pairs and squaring the 
result: 



a(e + e- -> NN) oc 



qeN 



(2) 



where (s) denotes the amplitude at E^ M = s for making the baryon 
N with a given primary flavor q, which is determined by the baryon wave 
functions. 

Since the wave functions of the baryon octet have a mixed symmetry, 
the amplitudes (s) tend to be highly asymmetric in specific models. For 
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Figure 1: Comparison of the cross sections for e + e~~ — > nn and pp in the 
threshold region E C m ~ 2 GeV. In the case of e + e~ — > pp, the direct- channel 
data are combined with the data for the time-reversed reaction pp — > e + e~ 
(marked by x ). The dash-dotted and dotted lines denote the average and l-o 
error bars, respectively, for the pp and nn data sets. 
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example, in the Chernyak-Zhitnitsky proton wave function ||, the u quark 
dominates, i.e., a v u = 0(1), a p d <C 1 and similarly 0% = 0(1), a™ <C 1. In such 
a limiting case we have 

a{e + e~ -> gp) Qj = 

a(e+e- -»■ nn) <2d ' 

While this is an extreme case, on general grounds we expect that the u 
contribution dominates in the proton and the d in the neutron, so that 
cr(e + e~ — > pp)/<j(e + e~ — > nn) ^> 1 at large momentum transfers. 

We find it puzzling that the experimental ratio ([!]) is apparently below 
unity when Eq M = s ~ 4 GeV 2 , whereas the ratio should be much larger than 
unity at both larger (timelike) and smaller (spacelike) momentum transfers. 
Clearly, the mechanism at work here is qualitatively different from those 
responsible for the above intuition. 

The lack of a conventional theoretical explanation is part of the moti- 
vation for the proposed new asymmetrical e~e + high-statistics collider at 
SLAC for the regime 1.4 < y/s < 2.5 GeV |J. This machine will yield 
high-precision data on baryon production in e~e + annihilation at thresh- 
old, providing a check on the FENICE data and an accurate benchmark for 
testing possible theoretical explanations. 

The first thing one must realize is that even though q 2 ^ 4m^ ^> Aq CD , 
the process is highly nonperturbative. This is because the 'extra' kinetic 
energy available to the quarks is very small. 

Our approach to this puzzle is based on thinking about the time- 
reversed processes: NN — > e + e~. These may be viewed as two-step pro- 
cesses, with a coherent meson state serving as an intermediate state. One pos- 
sible motivation for this picture might be provided by the Skyrme model || 
H, according to which baryons appear as solitons in a purely bosonic chiral 
Lagrangian. This model is formally justified as a low-energy approximation 
to large- A^ c QCD [|K| p| , and is known to provide a good description of 
many low-energy properties of baryons: see [T2|, [HJ for reviews. Skyrmion- 



anti-Skyrmion annihilation provides fl4]|-|17| a fairly accurate description 
of low-energy baryon-antibaryon annihilation. Just after the Skyrmion and 
anti-Skyrmion touch, they 'unravel' each other, and a coherent classical pion 
wave emerges as a burst that takes away energy and baryon number as quickly 
as causality permits. A specific parametrization of the initial pion configu- 
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ration is [JIB] : 

F(r,t = 0) =h^— 2 e- r ' a } (4) 

where F is the profile of the chiral field, U = exp[i r • f F(r, t) ], a is a range 
parameter, /i is chosen so that the total energy is that of the NN pair, and 
the form of F guarantees that the pion configuration has zero net baryon 
number. This crude model has been shown |TJ| to reproduce satisfactorily 
the inclusive single-pion spectrum in pp annihilation at rest and the branching 
ratios for multi-pion final states. 

The details of this specific configuration are unimportant for our pur- 
poses: what is important is that the data are not inconsistent with such a 
model. Indeed, although the Skyrme model provides some motivation for our 
approach, it is not even essential for our purpose. What is important is that a 
single intermediate state should dominate the two-step NN — > e + e~ process. 
This could, for example, equally well be a single intermediate J PC = 1 
resonant meson state. 

To be more precise, since NN annihilation is a strong-interaction process, 
one must consider separately the 1 = 1 and 1 = channels. Accurately 
stated, our key assumption is that both of these channels are dominated by 
single states. These might be some excited p* and uj* mesons, for example, 
just as well as coherent pion configurations (||) with 1 = 1 and 1 = 0. 

With this picture in mind, we write the I = 1, NN — > e + e~ annihilation 
amplitudes as A%, e^Ao, where the overall phase is irrelevant, A\ and Aq are 
relatively real, and a is the relative phase between the 1 = 1 and 1 = 
amplitudes. We then have 



./ 



<r(e + e — > pp) 



aye^e — > nn) 



A x + e ia A c 



Ai - e ia A Q 



(5) 



It is apparent from (^|) that cr(e + e~ — > nn)/o~(e + e~ — > pp) ~ 1 if either 
A\ 3> Aq or vice versa. 

Remarkably, there is evidence from both e + e~ and NN annihilations that 
1 = 1 final states dominate by large factors. 

The clearest evidence comes from e + e~ — > nrc, where it is found by mea- 
suring final states with even and odd numbers of pions respectively that 

a(e+e- - (2m)*) g ^ ^ ^ 2 ^ (g) 



cr(e+e~ -> (2m + l)n) 
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Figure 2: Cross sections for e + e~ — > multi-pion final states, for E C m ~ 
2 GeF Jl5|/ . VFe note i/ie dominance of I = 1 /ina/ states with even numbers 
ofpions by about an order of magnitude over I = states with an odd number 
of pions. 
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as seen in Fig. ^.[| At these energies, we expect most final states created 
by e + e~ — > ss to contain KK pairs, so that (^j) corresponds to the non-ss 
initial states we expect to dominate in NN annihilation. The value (^) is 
similar to that found at lower energies, where T(p — > e + e~) ~ 9 x V(u — > 
e + e~), in agreement with naive quark models. The fact that the ratio a (I = 
l)/a(I = 0) continues to be large at higher energies is consistent with ideas 
of generalized vector meson dominance. The data from NN annihilations are 
less clear. Theoretically, there are various calculations and other suggestions 
that the cross sections for NN annihilations into 1 = 1 and 1 = may 
be similar. Experimentally, several initial states contribute, including 1 So, 
3 S/D! and various P waves, but we are interested only in the J PC = 1 
3 S/Di initial states. It is in principle possible to distinguish different initial 
states by comparing annihilations in gas and liquid, as has been done in the 
analysis of OZI-violating final states, but we are unaware of a comparable 
analysis of multi-pion final states. Because the initial state is a mixture with 
different G = ±1, it is not possible to separate 1 = 1 from 1 = simply by 
counting pions, as was the case in e + e~ annihilation. 

The most convincing experimental information known to us comes from 
an analysis of NN — > KK. By comparing the rates for pp — > K + K~ and 
pp — > K°TC, it has been possible to extract 



A( 3 S/D 1 -> KK) I= 



A(*S/D 1 - KK) I=0 



5 to 10, 



(7) 



which is comparable to the corresponding ratio (^|) in e + e~ annihilation. The 
fact that 1=1 dominates over / = in the ratio (0) is consistent with the 
hypothesis that most of the KK final states are created by uu and dd pairs 
in the initial NN state, with a ss pair popping out of the vacuum. The ratio 
(0) would be small if primary ss pairs dominated. 

We now use the experimental information on the dominance by the 1 = 1 
channel in both e + e" (|6|) and 3 S/Di annihilation (0) in a quantitative analysis 
of the e + e~ — > NN production / ratio (Q). Defining e = A /Ai, we can 
rewrite (^]) as 

2 



/ 



1 + e ia e 



(8) 



^The five-pion final state is predominantly wirir. The cross section in Fig, p corresponds 
to the final state a>7r + 7r~; for the total ujittt one should multiply it by 1.5 |19|. 
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It is clear that the zero-momentum-transfer limit cr(e + e~ — > pp) ^> (e + e~ — > 
nn) is obtained in the limit e — > 1, a — > 0, and that the high- momentum- 
transfer limit cr(e + e~ — > pp) ~4x (e + e~ — > rm) is obtained in the limit e — > 
1/3, a — > 0. In order to estimate e = Ai/A , our assumption of dominance 
in each isospin channel by a single state (either a coherent multi-pion state 
(fD or a generalized vector meson V*) tells us that 



a{e+e ~ _> (/ = o)) 

x 



^iwv - (J = 0)) 



^ a{e+e~ -»■ (/ = 1)) \| a(iViV -»■ (/ = 1)) ' 
Inserting the experimental indications (|^,0) into (P), we estimate that 



The top end of this range seems to us quite conservative, whereas the lower 
end surely requires more justification from NN annihilation data. In the 
following numerical analysis, we keep e general, but focus extra attention on 
the limits e = 1/3 and 1/10. 

It is apparent from @ that / ~ 1 is possible for any value of e, for a 
restricted range of the relative phase a ~ n/2. However, the allowed range 
of a is extended if e is small. It is easy to see that / lies in a narrow range 
Af around unity if a falls within the following range: 

Act ~ -J-. (11) 
4e 



It is apparent that / ~ 1 for all values of a if e is small, as suggested fllQ|) by 
the available data on e + e~ and NN annihilation. 

The quantitative behaviour of / as a function of0<e<l/2 and —n < 
a < it is shown in Fig. |[ Displayed explicitly is the region of the (e, a) plane 
where / falls within the experimental range (|I]). We see that this range 
favours |a| > ir/2, whatever the value of e. Fig. f| displays projections of 
Fig. ^ for the two limiting values e = 1/3 and 1/10. The allowed range (j^) 
of / and the corresponding ranges of a are also shown. 

We conclude that the a priori puzzling large experimental value of the 
ratio cx(e + e _ — > nn)/a(e + e~ — > pp) can be understood qualitatively. This is 
relatively easy if the 1 = 1 amplitude dominates over the 1 = amplitude, 
as suggested by the available data on e + e~ and NN annihilation and our 
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Figure 3: Three-dimensional plot of the cross-section ratio f = er(e + e~ • 
pp)/a(e + e~ — > fin) as a function of e and a, indicating the region where 
falls within the range 
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Figure 4: Two-dimensional plot of the cross-section ratio f = a(e + e~ — > 
pp)/o-(e + e~ — > nn) as a function of a for e = 1/3 (dot-dash curve) and 
e = 1/10 (continuos curve), indicating the ranges where f falls within fifty. 
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assumption of dominance by a single coherent state in each isospin channel. 
The specific range ([!]) can be understood quantitatively if the 1 = 1 and 
1 = amplitude have a large relative phase a. We are not in a position to 
judge the plausibility of such a large value of a, from either an experimental 
or a theoretical point of view. It would be interesting to make tests of this 
possibility. 

We comment finally on the surprisingly large value of the ratio cr(77 — » 
AA)/cx(77 -> pp) observed by the CLEO Collaboration |g (see also 
for related experimental work). CLEO find that (7(77 — > pp) ~ (1(77 — > AA) 
close to threshold, which seems analogous to FENICE result for the nn/pp 
ratio. For a given quark flavor, the perturbative amplitude for baryon- 
antibaryon production in the photon-photon reaction scales like the quark 
charge squared, compared with the linear dependence of the amplitudes on 
the quark charge in the e + e~ case discussed earlier. Thus one might naively 
expect the ratio cr(AA) /a(pp) to be even smaller than the corresponding 
perturbative prediction for a(nn)/a(pp) in e + e~. It would be interesting to 
approach this puzzle from a point of view similar to that adopted in this 
paper. However, the situation in 77 collisions is more complicated, because 
of the wider range of possible spin and isospin states. Also, the informa- 
tion available on the isospin and spin decomposition is sparse compared with 
that in e + e~ annihilation, which we used above. Data for 77 — > nn close to 
threshold might cast light on the 17(77 ~^ AA)/cr(77 — > pp) puzzle, but are 
not yet available. 

We have proposed in this paper a simple model that is able to accommo- 
date the suprisingly large observed value of the ratio <r(e + e" — > nn) / <j(e + e~ — 
pp) . Our suggestion is based on a simple two-step approach, in which a single 
intermediate state with 1 = 1 dominates over 1 = 0. This dominant inter- 
mediate state could be motivated by a Skyrmion-anti-Skyrmion picture, or 
could be some excited p* resonance. Our model could be tested by further 
measurements of the ratios of different isospin amplitudes in e + e~ and NN 
annihilation, and suggests a relatively large phase difference between 1=1 
and 1 = amplitudes. We look forward to more experimental data bearing 
on these issues, for example from a new low-energy e + e~ collider ||. 
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